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We develop the hydrodynamic theory of collinear spin currents coupled to magnetization dynamics 
in metallic ferromagnets. The collective spin density couples to the spin current through a U(l) 
Berry-phase gauge field determined by the local texture and dynamics of the magnetization. We 
determine phenomenologically the dissipative corrections to the equation of motion for the electronic 
current, which consist of a dissipative spin-motive force generated by magnetization dynamics and 
a magnetic texture-dependent resistivity tensor. The reciprocal dissipative, adiabatic spin torque 
on the magnetic texture follows from the Onsager principle. We investigate the effects of thermal 
fluctuations and find that electronic dynamics contribute to a nonlocal Gilbert damping tensor in 
the Landau-Lifshitz-Gilbert equation for the magnetization. Several simple examples, including 
magnetic vortices, helices, and spirals, are analyzed in detail to demonstrate general principles. 
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I. INTRODUCTION 

The interaction of electrical currents with magnetic 
spin texture in conducting ferromagnets is presently a 
subject of active research. Topics of interest include 
current-driven magnetic dynamics of solitons such as do- 
main walls and magnetic vortices EEEE as we u ^ t ne 
reciproca l process of v oltage generation by magnetic 
dynamics! 5 | 6 | 7 | 8 | 9 | 1 ° | 11 1 12 1 This line of research has been 
fueled in part by its potential for practical applications 
to magnetic memory and data storage devicesPSl Funda- 
mental theoretica l inter est in the subject dates back at 
least two decades? 5 * 6 * 14 ! It was recognized early orP that 
in the adiabatic limit for spin dynamics, the conduction 
electrons interact with the magnetic spin texture via an 
effective spin-dependent U(l) gauge field that is a local 
function of the magnetic configuration. This gauge field, 
on the one hand, gives rise to a Lorentz force due to 
"fictitious" electric and magnetic fields and, on the other 
hand, mediates the so-called spin-transfer torque exerted 
by the conduction electrons on the collective magnetiza- 
tion. An alternative and equivalent view is to consider 
this force as the result of the Berry phase^ accumulated 
by an electron as it propagates through the ferromagnct 
with its sp in aligned with the ferromagnetic exchange 
field ISliSGSI j n the standard phenomenological formalism 
based on the Landau-Lifshitz-Gilbert (LLG) equation, 
the low-energy, long-wavelength magnetization dynamics 
are described by collective spin precession in the effective 
magnetic field, which is coupled to electrical currents via 
the spin-transfer torques. In the following, we develop 
a closed set of nonlinear classical equations governing 
current-magnetization dynamics, much like classical elec- 
trodynamics, with the LLG equation for the spin-texture 
"field" in lieu of the Maxwell equations for the electro- 
magnetic field. 

This clcctrodynamic analogy readily explains various 
interesting magnetoclcctric phenomena observed recently 
in ferromagnetic metals. Adiabatic charge pumping by 



magnetic dynamics^ can be understood as the gener- 
ation of electrical currents due to the fictitious electric 
field™ In addition, magnetic textures with nont rivial 
topology exhibit the so-called topological Hall effect J2ES1 
in which the fictitious magnetic field causes a classical 
Hall effect. In contrast to the classical magnetoresis- 
tance, the flux of the fictitious magnetic field is a topo- 
logical invariant of the magnetic texture P 

Dissipative processes in current-magnetization dynam- 
ics are relatively poorly understood and are of central 
interest in our theory. Electrical resistivity due to quasi- 
one-dimensional (ID) domain w alls and spin spirals have 
been calculated microscopically.^ 1 * 22 * 23 * More recently, a 
viscous coupling between current and magnetic dynam- 
ics which determines the strength of a dissipative spin 
torque in the LLG equation as well the reciprocal dis- 
sipative spin electromotive force generated by magnetic 
dynamics, called the "/3 coeffici ent,''^ was also calcu- 
lated in microscopic approaches! 3 * 24 * 25 -^ Generally, such 
first-principles calculations are technically difficult and 
restricted to simple models. On the other hand, the num- 
ber of different forms of the dissipative interactions in the 
hydrodynamic limit arc in general constrained by sym- 
metries and the fundamental principles of thermodynam- 
ics, and may readily be determined phenomenologically 
in a gradient expansion. Furthermore, classical thermal 
fluctuations may be easily incorporated in the theoretical 
framework of quasistationary nonequilibrium thermody- 
namics. 

The principal goal of this paper is to develop a (scmi- 
phenomcnological) hydrodynamic description of the dis- 
sipative processes in electric flows coupled to magnetic 
spin texture and dynamics. In Ref. [TTJ we drew the anal- 
ogy between the interaction of electric flows with quasis- 
tationary magnetization dynamics with the classical the- 
ory of magnetohydrodynamics. In our "spin magnetohy- 
drodynamics," the spin of the itinerant electrons, whose 
flows are described hydrodynamically, couples to the lo- 
cal magnetization direction, which constitutes the col- 
lective spin-coherent degree of freedom of the electronic 
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fluid. In particular, the dissipative (3 coupling between 
the collective spin dynamics and the itinerant electrons 
is loosely akin to the Landau damping, capturing cer- 
tain kinematic equilibration of the relative motion be- 
tween spin-texture dynamics and electronic flows. In our 
previous paper | we considered a special case of incom- 
pressible flows in a ID ring to demonstrate the essential 
physics. In this paper, we establish a general coarse- 
grained hydrodynamic description of the interaction be- 
tween the electric flows and textured magnetization in 
three dimensions, treating the itinerant electron's degrees 
of freedom in a two-component fluid model (correspond- 
ing to the two spin projections of spin-1 /2 electrons along 
the local collective magnetic order) . Our phenomenology 
encompasses all the aforementioned magnetoelectric phe- 
nomena. 

The paper is organized as follows. In Sec. [IT] we use a 
Lagrangian approach to derive the semiclassical equation 
of motion for itinerant electrons in the adiabatic approx- 
imation for spin dynamics. In Sec. |III[ we derive the 
basic conservation laws, including the Landau-Lifshitz 
equation for the magnetization, by coarse-graining the 
single-particle equation of motion and the Hamiltonian. 
In Sec. |IV[ we phenomenologically construct dissipative 
couplings, making use of the Onsager reciprocity princi- 
ple, and calculate the net dissipation power. In particu- 
lar, we develop an analog of the Navier-Stokes equation 
for the electronic fluid, focusing on texture-dependent 
effects, by making a systematic expansion in nonequi- 
librium current and magnetization consistent with sym- 
metry requirements. In Sec. |Vj we include the effects of 
classical thermal fluctuations by adding Langevin sources 
to the hydrodynamic equations, and arrive at the central 
result of this paper: A set of coupled stochastic differ- 
ential equations for the electronic density, current, and 
magnetization, and the associated white-noise correlators 
of thermal noise. In Sec. VI we apply our results to 



given by (see appendix |A| for details): 



special examples of rotating and spinning magnetic tex- 
tures, calculating magnetic texture resistivity and mag- 
netic dynamics-generated currents for a magnetic spiral 



and a vortex. The paper is summarized in Sec. VII and 
some additional technical details, including a microscopic 
foundation for our semiclassical theory, are presented in 
the appendices. 



II. QUASIPARTICLE ACTION 

In a ferromagnet, the magnetization is a symmetry- 
breaking collective dynamical variable that couples to the 
itinerant electrons through the exchange interaction. Be- 
fore developing a general phcnomcnological framework, 
we start with a simple microscopic model with Stoner in- 
stability, which will guide us to explicitly construct some 
of the key magnetohydrodynamic ingredients. Within a 
low-temperature mean-field description of short-ranged 
electron-electron interactions, the electronic action is 



S = / dtd 3 rif; 



ihdt 



2m P 



-V 2 - 
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Here, A(r, i) is the ferromagnetic exchange splitting, 
m(r, t) is the direction of the dynamical order param- 
eter defined by h(^ orxjj) /2 = p s m, p s is the local spin 
density, and ip(r,t) is the spinor electron field operator. 
For the short-range repulsion U > discussed in ap- 
pendix [Aj A(r,t) = 2Up s (r,t)/h and 0(r,t) = Up(r,t), 
where p — (ip^ip) is the local particle number density. 
For electrons, the magnetization M is in the opposite di- 
rection of the spin density: M = 7p s m, where 7 < is 
the gyromagnetic ratio. Close to a local equilibrium, the 
magnetic order parameter describes a ground state con- 
sisting of two spin bands filled up to the spin-dependent 
Fermi surfaces, with the spin orientation defined by m. 
We will focus on soft magnetic modes well below the 
Curie temperature, where only the direction of the mag- 
netization and spin density are varied, while the fluctu- 
ations of the magnitudes are not significant. The spin 
density is given by p s = h{p + — p~)/2 and particle den- 
sity by p = P+ + P-, where p± are the local spin-up/down 
particle densities along m. p s can be essentially constant 
in the limit of low spin susceptibility. 

Starting with a nonrelativistic many-body Hamilto- 
nian, the action is obtained in a spin-rot ationally 
invariant form. However, this symmetry is broken by 
spin-orbit interactions, whose role we will take into ac- 
count phenomenologically in the following. When the 
length scale on which m(r, t) varies is much greater than 
the ferromagnetic coherence length l c ~ hvp/A, where 
vf is the Fermi velocity, the relevant physics is captured 
by the adiabatic approximation. In this limit, we start 
by neglecting transitions between the spin bands, treat- 
ing the electron's spin projection on the magnetization 
as a good quantum number. (This approximation will 
be relaxed later, in the presence of microscopic spin- 
orbit or magnetic disorder.) We then have two effec- 
tively distinct species of particles described by a spinor 
wave function ip' , which is defined by tp — U(lZ)?p' . Here, 
U(1Z) is an SU(2) matrix corresponding to the local spa- 
tial rotation lZ(r,t) that brings the z-axis to point along 
the magnetization direction: H(r,t)z — m(r, t), so that 
W (or ■ m)U = a z . The projected action then becomes: 



dt / d 3 r^ n 



(ihd t + a) 



(p A „ 

1 <j ? 

2 2 2 



{-ihV - a) 2 

2m e 

- f dtF[m] , (2) 



where 



F[m] = ^ J d 3 r(d 4 m) 2 



(3) 



is the spin-texture exchange energy (implicitly summing 
over the repeated spatial index i), which comes from the 
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terms quadratic in the gauge fields that survive the pro- 
jection. In the mean-field Stoner model, the ferromag- 
netic exchange stiffness is A — h 2 p/Am e . To broaden our 
scope, we will treat it as a phcnomenological constant, 
which, for simplicity, is determined by the mean electron 
density.^ The spin-projected "fictitious" gauge fields are 
given by 



a a (r,t) = ih(a\U^ d t U\a) 
a CT (r,i) = ih{a\U ] VU\a) 



(4) 



Choosing the rotation matrices U(m) to depend only on 
the local magnetic configuration, it follows from their 
definition that spin-er gauge potentials have the form: 

fl^-ftm-anm), a m = -9,m • , (5) 

where a" lon (m) = — ih(a\W d m U\a) . We show in Ap- 
pendix [B] the well known result (see, e.g., Ref. [27jl that 
a™ is the vector potential (in an arbitrary gauge) of 
a magnetic monopolc in the parameter space defined by 
m: 



9 m xar» = 9 ff m, 



(6) 



where q a = ah/2 is the monopolc charge (which is ap- 
propriately quantized) . 

By noting that the action Q is formally identical to 
charged particles in electromagnetic field, we can imme- 
diately write down the following classical single-particle 
Lagrangian for the interaction between the spin-er elec- 
trons and the collective spin texture: 

L a (r, r, t) = — | h r ■ a ff (r, t) + a a (r, t) , (7) 

where r is the spin-er electron (wave-packet) velocity. To 
simplify our discussion, we are omitting here the spin- 
dependent forces due to the self-consistent fields <fi(r,t) 
and A(r,t), which will be easily reinserted at a later 
stage. See Eq. (§91. 



The Euler-Lagrange equation of motion for v = 
r derived from the single-particle Lagrangian ([7]), 
(d/dt)(dL a /dr) = dL a /dr, gives 



m e v = q a (e + v x b) 



(8) 



The fictitious electromagnetic fields that determine the 
Lorcntz force are 

q a ei = dida - d t a ai = q a m ■ (d t m x d l m) , 



q<rh = e'^djdak = q a —m ■ (<9 fc m x djva) . (9) 

They are conveniently expressed in terms of the tensor 
field strength 

Q<rUv = d^a^ - d v a ail = q a m ■ (d v m x d^m) (10) 

by e« = fio an d bi = e 1 ^ fjk/2. e 1 ^ is the antisymmet- 
ric Levi-Civita tensor and we used four-vector notation, 



ijk 



defining <9 M = (9 t ,V) and a afl = (a a ,a a ). Here and 
henceforth the convention is to use Latin indices to de- 
note spatial coordinates and Greek for space-time coor- 
dinates. Repeated Latin indices k are, furthermore, 
always implicitly summed over. 



III. SYMMETRIES AND CONSERVATION 
LAWS 

A. Gauge invar iance 

The Lagrangian describing coupled electron transport 
and collective spin-texture dynamics (disregarding for 
simplicity the ordinary electromagnetic fields) is 



L(r p ,v p ;m,9 M m) 



E 

p 

E 



m e v; 



A 



v„ ■ a„ + a„ 



d 3 r{d l xnf 



^Xm) ~\j d 3 r(d t m) 2 . (11) 



wp 1 = (l,v p ), v p = r, and cr here is the spin of indi- 
vidual particles labelled by p. The resulting equations 
of motion satisfy certain basic conservation laws, due to 
spin-dependent gauge freedom, space-time homogeneity, 
and spin isotropicity. 

First, let us establish gauge invariance due to an ambi- 
guity in the choice of the spinor rotations U(r,t) — > UU' . 
Our formulation should be invariant under arbitrary di- 
agonal transformations U' = e~ l f and U' — e^ t9,J '' 2 on 
the rotated fermionic field ip', corresponding to gauge 
transformations of the spin-projected theory: 



<5a CTAl = hdfif and 5a all = , 



(12) 



respectively. The change in the Lagrangian density is 
given by 



5C = j> l d ll f and 8C = tfd ll g, 



(13) 



respectively, where j = j + + j_ and j s = h(j + — j-)/2 
are the corresponding charge and spin gauge currents. 
The action S = j dtd 3 rC is gauge invariant, up to sur- 
face terms that do not affect the equations of motion, 
provided that the four-divergence of the currents vanish, 
which is the conservation of particle number and spin 
density: 



/?+V-j = 0, p s + V-j s =0. 



(14) 



(The second of these conservation laws will be relaxed 
later.) Here, the number and spin densities along with 
the associated flux densities are 



p = E n p = p+ + p- 

p 

j = n pV p = pv , 



(15) 
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and 



Pa = ^2q<rn p = ~{p+ - p-) 
P 

is = ^qaUpVp EE p s V s , 



(16) 



where n p = S(r — r p ) and u v = ± for spins up and down. 
In the hydrodynamic limit, the above equations deter- 
mine the average particle velocity v and spin velocity 
v s , which allows us to define four- vectors — (p,pv) 
and j£ — (p s ,jO s v s ). Microscopically, the local spin- 
dependent currents are defined, in the presence of electro- 
magnetic vector potential a and fictitious vector potential 
a CT , by 



where e < is the electron charge. 



B. Angular and linear momenta 



(17) 



where the angular-momentum stress tensor is defined by 

Ily = p s v s] mi - A(m x . (23) 

Notice that this includes both quasiparticle and collective 
contributions, which stem respectively from the trans- 
port and equilibrium spin currents. 

The Lorentz force equation for the electrons, Eq. pj, 
in turn, leads to a continuity equation for the kinetic 
momentum densityP To see this, let us start with the 
microscopic perspective: 



(24) 



Using the Lorentz force equation for the second term, we 
have: 

m e ^ n pVp = Y 1<r n p( e i + eVkb kV p] ) = ^ <la n pfiH V p 
P P P 

= p s m ■ (d t m x dim) + p s v S jm • (djm x c^m) 
= (dim) ■ (S m F) = -A(d % m) ■ (3 2 m) , (25) 



Our Lagrangian ( |11[ ) contains the dynamics of m(r) 
that is coupled to the current. In this regard, we note 
that the time component of the fictitious gauge poten- 
tial (B4|, a a — —Hdt(p(l — <J cos 9) /2, is a Wess-Zumino 
action that governs the spin-texture dynamics ! 4 * 6 * 28 ! The 
variational equation m x 5 m L = gives: 



p s (d t + v s • V)m + m x <5 m F = . 



To derive this equation, we used the spin-density con- 
tinuity equation (14) and a gauge-independent identity 



satisfied by the fictitious potentials: their variations with 
respect to m are given by 



where 



S m a afl (m, d^m) = q a m x c^m . 



dm ^"^d(d ll m) 



(19) 
(20) 



One recognizes that Eq. (18 1 is the Landau-Lifshitz (LL) 



equation, in which the spin density precesses about the 
effective field given explicitly by 



h == 5 m F = -Adfm . 



(21) 



Equation ( |18[ ) also includes the well-known reactive spin 
torque: r = (j s • V)mP which is evidently the change 
in the local spin-density vector due to the spin angular 
momentum carried by the itinerant electrons. One can 
formally absorb this spin torque by defining an advective 
time derivative D t = dt + v s • V, with respect to the 
average spin drift velocity v s . 

Equation ( 18 I may be written in a form that explicitly 



expresses the conservation of angular momentum! 27 ! 29 -) 



utilizing Eq. ([18]) to obtain the last line. Coarse-graining 
the first term of Eq. p4|, in turn, we find: 



^p v p = ~ d J Y s ( r - v p) 



VrtiV'n 



~9j Y PaVaiVaj ■ 

(26) 

Putting Eqs. (|25| and (|26| together, we can finally write 



(18) Eq. (24) in the form 



rn 



■dt(pvi) + dj ^Tij + m e Y P<y 



where 



Tij — A 



(dim) ■ (djm) 



f(d k mf 



, (27) 



(28) 



is the magnetization stress tensorP 

A spin-dependent chemical potential p = K~ x p gov- 
erned by local density and short-ranged interactions can 
be trivially incorporated by redefining the stress tensor 
as 



T 



(29) 



In our notation, p = (p+, P-) T , p = (p+,p-) T and K is 
a symmetric 2x2 compressibility matrix in spin space, 
which includes the degeneracy pressure as well as self- 
consistent exchange and Hartree interactions. In general, 
Eq. ( 29 1 is valid only for sufficiently small deviations from 



the equilibrium density. 



Using the continuity equations ( 14 1 , we can combine 



the last term of Eq. ( 27 1 with the momentum density 
rate of change: 



dt(p s rm) + djUij 



0. 



(22) 



dt(p<jVai) + dj(p tT V ai V a j) = Pcr(dt + V CT • V)^ , (30) 
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which casts the momentum density continuity equation 
in the Euler equation form: 



0. 



(31) 



We do not expect such advective corrections to dt to 
play an important role in electronic systems, however. 
This is in contrast to the advective-likc time derivative 



in Eq. (18), which is first order in velocity field and is 



crucial for capturing spin-torque physics. 



C. Hydrodynamic free energy 

We will now turn to the Hamiltonian formulation and 
construct the free energy for our magnetohydrodynamic 
variables. This will subsequently allow us to develop a 
noncquilibrium thermodynamic description. The canon- 
ical momenta following from the Lagrangian (11) are 



TV 



dL 
dv p 
dL 
dm 



p 



da a 
' 3m 



p 



(32) 



Notice that for our translationally-invariant system, the 
total linear momentum 



p 



J d 3 r(Tt • V) 



m = m, 



(33) 



where we have used Eq. ^ to obtain the second equality, 
coincides with the kinetic momentum (mass current) of 
the electrons. The latter, in turn, is equivalent to the lin- 
ear momentum of the original problem of interacting non- 
relativistic electrons, in the absence of any real or ficti- 
tious gauge fields. See appendix [A] While P is conserved 
(as discussed in the previous section and also follows now 
from the general principles), the canonical momenta of 
the electrons and the spin-texture field, Eqs. (32 1, are 



not conserved separately. As was pointed out by Volovik 
in Ref. this explains anomalous properties of the lin- 
ear momentum associated with the Wess-Zumino action 
of the spin-texture field: This momentum has neither 
spin-rotational nor gauge invariance. The reason is that 
the spin-texture dynamics define only one piece of the 
total momentum, which is associated with the coherent 
degrees of freedom. Including also the contribution as- 
sociated with the incoherent (quasiparticle) background 
restores the proper gauge-invariant momentum, P, which 
corresponds to the generator of the global translation in 
the microscopic many-body description. 

Performing a Legendre transformation to Hamiltonian 



as a function of momenta, we find 



H [r p , p p ; m, ir] = v p ■ p p + J <fVm ■ n - L 

+ | J d 3 r(^m) 2 



E 



V 

E + F. 



(34) 



where E is the kinetic energy of electrons and F is the 
exchange energy of the magnetic order. As could be 
expected, E is the familiar single-particle Hamiltonian 
coupled to an external vector potential. According to 
a Hamilton's equation, the velocity is conjugate to the 
canonical momentum: v p = dH/dp p . We note that ex- 
plicit dependence on the spin-texture dynamics dropped 
out because of the special property of the gauge fields: 
rh • 9ma<T = a a . Furthermore, according to Eq. (19 1, we 
have m x 8 m E = (j s • V)m, so the LL Eq. (f l8| can be 
written in terms of the Hamiltonian ( 34 1 aa^ 



p s m + m x S m H = . 



(35) 



So far, we have included in the spin-texture equa- 
tion only the piece coupled to the itinerant electron de- 
grees of freedom. The purely magnetic part is tedious 
to derive directly and we will include it in the usual LL 
phenomenology. 29 To this end, we redefine 



F[m(r)] -► F + F' , 



(36) 



by adding an additional magnetic free energy F'[m(r)], 
which accounts for magnetostatic interactions, crystalline 
anisotropies, coupling to external fields, as well as energy 
associated with localized d or / orbitalsP^ Then the to- 
tal free energy (Hamiltonian) is H = E + F, and we in 
general define the effective magnetic field as the thermo- 
dynamic conjugate of m: h = S m H . The LL equation 
then becomes 



g s m + m x h = , 



(37) 



where g s is the total effective spin density. To enlarge 
the scope of our phenomenology, we allow the possibility 
that g s =/= p s . For example, in the s — d model, an extra 
spin density comes from the localized <i-orbital electrons. 
Microscopically, g s d t m term in the equation of motion 
stems from the Wess-Zumino action generically associ- 
ated with the total spin density. 

In the following, it may sometimes be useful to separate 
out the current-dependent part of the effective field, and 
write the purely magnetic part as h m = S m F, so that 



h = h m - m x (j s • V)m 



and Eq. (37 1 becomes 



g s m + (j s ■ V)m + m x h„ 



0. 



(38) 



(39) 
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For completeness, let is also write the equation of motion 
for the spin-er acceleration: 

rn e (d t + v CT • V)v al = q a [m ■ (d t m x dim) 

+ v a jin ■ (djin x dim)] — V/i CT , (40) 

retaining for the moment the advective correction to 
the time derivative on the left-hand side and reinserting 
the force due to the spin-dependent chemical potential, 
fx = K~ x p. These equations constitute the coupled re- 
active equations for our magneto-electric system. The 
Hamiltonian (free energy) in terms of the collective vari- 
ables is (including the elastic compression piece) 



H \p a , p CT ; m] = ^ J d 3r P- 



(p<t - a a ) 2 
2m P 



d s rp 2 K^p + F[m] , (41) 



where p<j = m e v ff + a ff is the spin-dependent momentum 
that is locally averaged over individual particles. 



D. Conservation of energy 

So far, our hydrodynamic equations are reactive, so 
that the energy (41 1 must be conserved: P = H = E + 



F = 0. The time derivative of the electronic energy E is 
m e p a v G v a + p a 



m e p a v a jV a j - dj{p a v (yj ) 



m e v% 



= / d 3 r ^ PoV<j 3 [m e (d t + ■ V) v aj + djp a ] 

J a 

= d 3 r^2 <l<?PaV<j • (e + v ff x b) 

= J d3? ' laPa^a ■ e = J cf 3 rj s ■ e . (42) 

The change in the spin-texture energy is given, according 
to Eq. pi, by 



F = J d 3 rm ■ S m F = J d 3 rm ■ h m 

= J d 3 rm ■ [g s m x rh + m x (j s • V)m)] 

= - / d 3 rj s ■ e . (43) 



The total energy is thus evidently conserved, P = 0. 
When we calculate dissipation in the rest of the paper, 
we will omit these terms which cancel each other. The 
total energy flux density is evidently given by 



m e v„ 



IV. DISSIPATION 

Having derived from first principles the reactive cou- 
plings in our magneto-electric system, summed up in 
Eqs. ( 39 )-( 41 1 , we will proceed to include the dissipa- 



tive effects phcnomcnologically. Let us focus on the lin- 
earized limit of small deviations from equilibrium (which 
may be spin textured), so that the advective correction 
to the time derivative in the Euler Eq. ( |40| , which is 
quadratic in the velocity field, can be omitted. To elimi- 
nate the quasiparticle spin degree of freedom, let us, fur- 
thermore, treat halfmetallic ferromagnets, so that p = p + 
and p s — qp, where q = h/2 is the electron's spinP^From 



Eq. (40 1, the equation of motion for the local (averaged) 
canonical momentum isP^l 



p = -jxb-Vjii, 
P 



(45) 



in a gauge where a a = 0, so that p = m e v — geP2 
p = p/K. The Lorentz force due to the applied (real) 
electromagnetic fields can be added in the obvious way 
to the right-hand side of Eq. ( |45| . Note that since we 
are now interested in linearized equations close to equi- 
librium, p in Eq. (45) can be approximated by its (ho- 
mogeneous) equilibrium value. 

Introducing relaxation through a phcnomcnological 
damping constant (Drude resistivity) 
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pr 



(46) 



where r is the collision time, expressing the fictitious 
magnetic field in terms of the spin texture, Eq. (45 1 be- 
comes: 



Pi = --(m x ^111) • (j • V)m - dip - jji . (47) 
P 



Adding the phcnomcnological Gilbert damping 3 to 



the magnetic Eq. (37 1 gives the Landau-Lifshitz-Gilbert 
equation: 



g s (m + am x m) = h x m , 



(48) 



where a is the damping constant. Eqs. (47 1 and 
(48 1, along with the continuity equation, p — —V • j, 



are the near-equilibrium thermodynamic equations for 
(p, p,m) and their respective thermodynamic conjugates 
(p, j, h) = (S p H, S p H 7 S m H). This system of equations of 
motion may be written formally as 




= r[m(r)] 




(49) 



The matrix T depends on the equilibrium spin texture 
m(r). By the Onsager reciprocity principle, T^m] = 
(44) SiSjTji[— m], where Sj = ± if the ith variable is even 
(odd) under time reversal. 
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In the quasistationary description of a nonequilibrium 
thermodynamic system, the entropy S[p, p, m] is for- 
mally regarded as a functional of the instantaneous ther- 
modynamic variables, and the probability of a given con- 
figuration is proportional to e s ^ kB . If the heat conduc- 
tance is high and the temperature T is uniform and con- 
stant, the instantaneous rate of dissipation P = TS is 
given by the rate of change in the free energy, P = H = 
Jd 3 rV: 



into a tensor form: 



V = —pp — h • m — j • p = ag s m 2 



7J 



(50) 



where we used Eq. (47 1 and expressed the effective field 



h as a function of m by taking mx of Eq. (48 1 



h = g s m x m — ag s m . 



(51) 



Notice that the fictitious magnetic field b does not con- 
tribute to dissipation because it does not do work. 

So far, there is no dissipative coupling between the 
current and the spin-texture dynamics, and the macro- 
scopic equations obey the global time-reversal symme- 
try. However, we know that dissipative couplings ex- 
ists due to the misalignment of the electron's spin with 
the collective spin texture and spin-texture resistivity! 3 ^ 22 ! 
Following Ref. QTJ we add these well-known effects phe- 
nomenologically by making an expansion in the equations 
of motion to linear order in the nonequilibrium quanti- 
ties m and j. To limit the number of terms one can write 
down, we will only add terms that are spin-rotationally 
invariant and isotropic in real space (which disregards, 
in particular, such effects as the angular magnetoresis- 
tance and the anomalous Hall effect). To second order in 
the spatial gradients of m, there are only three dissipa- 
tive phenomenological terms with couplings ?y, rj , and (3 
consistent with the above requireme nts , which could be 
added to the right-hand side of Eq. (47lP3The momen- 
tum equation becomes: 



pi = (m x dim) ■ (j • V)m - d t p - jji 

P 

— r/(9fem) 2 ji — rj'diia ■ (j • V)m — q(3m • c^m . (52) 



7ij[m] =Sij [7 + T](d k m) 2 ] + rfditn ■ djm 
+ -m ■ (dim x djm) . 



(53) 



The last term due to fictitious magnetic field gives a Hall 
resistivity. Note that j[m] = 7 T [— m], consistent with 
the Onsager theorem. We can finally write Eq. (47 1 as: 



Pi = - Jij [ia]jj - difi - qf3m ■ dim . 



(54) 



As was shown in Ref. QT1 since the Onsager relations 
require that r[m] = T[— m] T within the current/spin- 
texture fields sector, there must be a counterpart to the 
(3 term above in the magnetic equation, which is the well- 
known dissipative "/3 spin torque:" 

g s (m + am x m) = h x m — q(3m x (j • V)m . (55) 

The total dissipation V is now given by 

V =ag s m 2 + 2q(3m ■ (j • V)m + [7 + ^(^m) 2 ] j 2 
+ r/[(j- V)m] 2 

OLQs 

ag s 



^ag s 



m 



(j • V)m 



+ [7 + 77(<9 fc m) 2 ] f 



[(j • V)m] 



(56) 



The second law of thermodynamics requires the total dis- 
sipation to be positive, which puts some constraints on 
the allowed values of the phenomenological parameters. 
We can easily notice, however, that the dissipation (56 1 
is guaranteed to be positive-definite if 



77 + 7]' > 



ag s 



(57) 



which may serve as an estimate for the spin-texture re- 
sistivity due to spin dephasing. This is consistent with 
the microscopic findings of Ref. l23l 



THERMAL NOISE 



It is known that the "/3 term" comes from a misalignment 
of the electron spin with the collective spin texture, and 
the associated dephasing. It is natural to expect thus 
that the dimensionless parameter (3 ~ H/t s A, where t s 
is a characteristic spin-dephasing timeP The "77 terms" 
evidently describe texture-dependent resistivity, which 
is anisotropic with respect to the gradients in the spin 
texture along the local current density. Such term are 
also naturally expected, in view of the well-known giant- 
magnetoresistance effect in which noncollinear magne- 
tization results in electrical resistance. The microscopic 
origin of this term is due to spin-texture misalignment, 
which modifies electron scattering. 

The total spin-texture-dependent resistivity can be put 



At finite temperature, thermal agitation causes fluc- 
tuations of the current and spin texture, which are cor- 
related due to their coupling. A complete description 
requires that we supplement the stochastic equations of 
motion with the correlators for these fluctuations. It 
is convenient to regard these fluctuations as being due 
to the stochastic Langevin "forces" (5/j,,Sj,Sh.) on the 
right-hand side of Eq. ( 49 ) . The complete set of finite- 



temperature hydrodynamic equations thus becomes: 

P= -v], 

p + q(3rhiVmi = -7[m]j - V/j , 
q s (1 + amx)m = h x m — qf]m x (j ■ V)m . (58) 
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where (fl, j, h) = (jit + <5/x, j + 5j, h + <5h). The simplest 
(while possibly not most realistic) case corresponds to 
a highly compressible fluid, such that K — ► oo. In this 
limit, /i = /o/AT — > and the last two equations com- 
pletely decouple from the first, continuity equation. In 
the remainder of this section, we will focus on this special 
case. The correlations of the stochastic fields are given 
by the symmetric part of the inverse matrix T = — r _1 j^l 
which is found by inverting Eq. ( |58[ ) (reduced now to a 
system of two equations) : 



Substituting this p in Eq. ( 59 1 , we can easily find a closed 



j = -7 1 (P + AVm,) , 

h = g s m x m — ag s m — q/3(j ■ V)m . 



(59) 



Writing formally these equations as (after substituting j 
from the first into the second equation) 



-T[m(r)] 



P 
m 



(60) 



we immediately read out for the matrix elements 
T(r,r') = T(r)£(r-r'): 

r iidA T ) =(7 _1 )«' , 
^ji,h tl (r) =ql3(7~ 1 )ikd k m i , , 
T 'v-:K( r ) = - 9/5(7^ 1 )fci9 fe m i / , 
Tfe ( ,h.,(r) =ag s S u > + g s e vl k m k 

- {qf5) 2 {d k m i ){r%k'{d k ,m il ) . (61) 

According to the fluctuation-dissipation theorem, we 
symmetrize T to obtain the classical Langevin 
correlators!^ 



{Sji{r,t)5j it {v',t'))/T=g iV , 
(5j i {r,t)5h i ,{r\t'))/T=qpg' ik d k m v , 
(5h i {r,t)5h i ,{r\t'))/T=aQ s 5 iV 

- (qf3) 2 gkk'(d k mi)(d k >mi>) , 

where T = 2k B TS{r - r')5(t - t') and 

g = [T 1 + (TYyi , g' = IT 1 - (ff]/2 



(62) 



(63) 



are, respectively, the symmetric and antisymmetric parts 
of the conductivity matrix 7 -1 . The short-ranged, 5- 
function character of the noise correlations in space stems 
from the assumption of high electronic compressibility. 
Contrast this to the results of Ref . [TTJ for incompressible 
hydrodynamics. A presence of long-ranged Coulombic 
interactions and plasma modes would also give rise to 
nonlocal correlations. These are absent in our treatment, 
which disregards ordinary electromagnetic phenomena. 

Focusing on the microwave frequencies u> characteris- 
tic of ferromagnetic dynamics, it is most interesting to 
consider the regime where u> <C t -1 . This means that 
we can employ the drift approximation for the first of 
Eqs. (l59l: 



stochastic equation for the spin-texture field: 

g s (l + amx)rh + m x rrh = (h m + Sh) x m , (65) 
where we have defined the "spin-torque tensor" 
r = q 2 (T 1 )kk> (m x d k m - (3d k m) 



(m x d k /m + /3d k /m) 



(66) 



The antisymmetric piece of this tensor modifies the effec- 
tive gyromagnetic ratio, while the more interesting sym- 
metric piece determines the additional nonlocal Gilbert 
damping: 



„ r + r T g 2 - 

a = — = — G, 

2p s Qs 



(67) 



where 



G =g kk ' [(m x d k m) ® (m x <9 fe 'm) - /3 2 d k m ® 9 fc /m] 
+ fig' kk : [(m x d k va) ® d k >m - d k m ® (m x d k >m)] 



In obtaining Eq. (|65|) from Eqs. (59 1, we have separated 



the reactive spin torque out of the effective field: h = 
h m -(jmx(j-V)m. (The remaining piece h m thus reflects 
the purely magnetic contribution to the effective field.) 
The total stochastic magnetic field entering Eq. (|65|, 



Sh = Sh + qm x (Sj ■ V)m , 



(69) 



captures both the usual magnetic Brown noise^ Sh 
and the Johnson noise spin-torque contribution's! Shj — 
qmx ((5j-V)m that arises due to the substitution j = j— <5j 
in the reactive spin torque q(j ■ V)m. Using correla- 
tors ( [62] ), it is easy to show that the total effective field 
fluctuations <5h are consistent with the nonlocal effec- 



tive Gilbert damping tensor ( 68 1 , in accordance with the 
fluctuation-dissipation theorem applied directly to the 
purely magnetic Eq. (651. 



To the leading, quadratic order in spin texture, we can 
replace g kk > — > S kk > /j and g' kk , — > in Eq. (68 1. This ad- 
ditional texture-dependent nonlocal damping (along with 
the associated magnetic noise) is a second-order effect, 
physically corresponding to the backaction of the magne- 
tization dynamics-driven current on the spin texture 
It should be noted that in writing the modified LLG 
equation ( 55 1 , we did not systematically expand it to 



include the most general phenomenological terms up to 
the second order in spin texture. We have only included 
extra spin-torque terms, which are required by the On- 



sager symmetry with Eq. ( 52 1 . The second-order Gilbert 



damping ( 68 1 was then obtained by solving Eqs. ( 52 ) and 
(551 simultaneously. (Cf. Refs. ITTI40I ') This means in 
particular, that this procedure does not capture second- 
order Gilbert damping effects whose physical origin is 
unrelated to the longitudinal spin-transfer torque physics 
studied here. One example of that is the transverse spin- 



P, 



m e Vi - qei 



-qei = qm ■ (m x dim) . (64) pumping induced damping discussed in Refs. ST] 
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VI. EXAMPLES 



and a [3 force 



A. Rigidly spinning texture 

To illustrate the rj resistivity terms in the electron's 
equation of motion (52 1, we first consider ID textures. 



Take, for example, the case of a ID spin helix m(z) 
along the z axis, whose spatial gradient profile is given by 
d z m = kz x m, where k is the wave vector of the spatial 
rotation and mlz. See Fig. [T] It gives anisotropic re- 



sistivity in the xy plane, 
Jv). 



, and along the z direction, 



rj(d z m.y 



TjK 



rfP = {v + t]')k 2 



(70) 




FIG. 1: (Color online) The transverse magnetic helix, d z m = 
kz x m, with texture-dependent anisotropic resistivity (70 1. 
We assume here translational invariance in the transverse (xy) 
directions. Spinning this helix about the vertical z axis gen- 
erates the dissipative electromotive forces which is spa- 
tially uniform and points everywhere along the z axis. A 
magnetic spiral, d z m — K(p x m = k6, spinning around the z 
axis, on the other hand, produces a purely reactive electromo- 
tive force e z , as discussed in the text, which is oscillatatory 
in space along the z axis. 

The fictitious electric field and dissipative j3 force re- 
quire magnetic dynamics. A general texture globally ro- 
tating clockwise in spin space in the xy plane according 
to m = —loz x m (which may be induced by applying a 
magnetic field along the z direction) generates an electric 
field 

e, = (m x m) • d{m — — to(m x z x m) • c^m 
= —bjdim z = —uidi cos 6 (71) 



= —/3m ■ dim — (3loz ■ (m x dim) 
= flusm 2 9diip , 



(72) 



where (9, tp) denote the position-dependent spherical an- 
gles parametrizing the spin texture. The reactive force 
( 71 1 has a simple interpretation of the gradient of the 
Berry-phase^ accumulation rate [which is locally deter- 
mined by the solid angle subtended by m(t)]. In the 
case of the transverse helix discussed above, 9 = tt/2, 
(p = kz — Lot, so that e 2 = while /i' 3 ' = —Pluk is finite. 

As an example of a dynamical texture that does not 
generate f ^ while producing a finite e, consider a spin 
spiral along the z axis, described by d z m = Ktpxm = k6, 
and rotating in time in the manner described above. It is 
clear geometrically that the change in the spin texture in 
time is in a direction orthogonal to its gradients in space. 



Specifically, 9 = kz, tp = —tot, so that /. 
electric field is oscillatory, e z = LOKsm.9 



while the 



B. Rotating spin textures 

We show here that a vortex rotating about its core in 
orbital space generates a current circulating around its 
core, as well as a current going radially with respect to 
the core. Consider a spin texture with a time depen- 
dence corresponding to the real-space rotation clockwise 
in the xy plane around the origin, such that m(r, t) — 
m(r(f), 0) with r = wzxr = ujr<fi, where we use polar co- 
ordinates (r, </>) on the plane normal to the z axis in real 
space [to be distinguished from the spherical coordinates 
(9, if) that parametrize m in spin space], we have 



rh = (f • V)m = Lod^m . 



(73) 



For m(r, <f>) in polar coordinates, the components of the 
electric field are, 



e r = cum ■ (d^m x d r m) , = , 
while the components of the (3 force are 

(d^m) 2 



03) 



-0U>- 



(74) 



(75) 



In order to find the fictitious electromagnetic fields, we 
need to calculate the following tensors (which depend on 
the instantaneous spin texture): 



bij = m • (dim x djm) = sin 9(di9djtp — dj9di<p) , 



dij = dim ■ djm 



di9dj9 + sin 2 9ditpdjip . 



(76) 



As an example, consider a vortex centered at the ori- 
gin in the xy plane with winding number 1 and positive 
polarity, as shown in Fig. [2] Its angular coordinates are 
given by 



tp = (0 + LOt) + 



2 ' 



9(r) 



(77) 
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where 4> = arg(r) and 9 is a rotationally invariant func- 
tion such that 9 — > as r — > and — > it/2 as r — > oo. 
Evaluating the tensors in equation ([76]) for this vortex in 
polar coordinates gives d rr = (d r 9y% — (sin#/r) 2 , 
d rt p = 0, and b r< p — ~(d r cos9)/r. The radial electric 
field is then given by 



-Ljrb r 



cod r cos 9 . 



The P force is in the azimuthal direction: 



-j3ujrd, 



4>4> 



sin 2 (9 



(78) 



(79) 



We can interpret this force as the spin texture "dragging" 
the current along its direction of motion. Notice that the 
forces in Eqs. (78 1 and (79) are the negative of those in 



Eqs. ( 71 ) and (72 1, as they should be for the present case, 



since the combination of orbital and spin rotations of our 
vortex around its core leaves it invariant, producing no 
forces. 



V 



rtttttt 

■s / / / ; rt 



FIG. 2: Positive-polarity magnetic vortex configuration pro- 
jected on the xy plane, m has a positive (out-of-plane) z 
component near the vortex core. Rotating this vortex about 
the origin in real space generates the current in the xy plane 
shown in Fig. [3] 



The total resistivity tensor (53 1 is (in the cylindrical 
coordinates) 



7 = 7 + n (d rr + d^) + rfd+ -b = ( _^ 7± 

p \ 7-1 I4> 

where 



(80) 



J r = J+(v + v')(d r 9) 2 + Ti 



sin 9 



l4> = 1 + i 1 {d r 9) 2 + {i 1 + i 1 ') (_j 

q d r cos 9 
p r 



7-L = 



(81) 



Here, the two diagonal components, 7 r and 7^, describe 
the (dissipative) anisotropic resistivity, while the off- 
diagonal component, captures what is called the 
topological Hall effect 

In the drift approximation, Eq. (|64|), the current- 
density field j = j r r - 



j = 1 1 q(e- 



qw7 



qui sin 9 



Eq. ([641 

- jgd is given by 

d r COS 9 

-/3sin 2 9/r 



lrl<t> 



7* 
7i 



-7-1 

7r 



d r 9 
(3sm9/r 



(82) 



More explicitly, we may consider a profile 9 = 7r(l — 
e~ r / a )/2, where a is the radius of the vortex core. The 
corresponding current (82 1 is sketched in Fig. [3] 
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FIG. 3: We plot here the current in Eq. (82 1 (all parameters 



set to 1). Near the core, the current spirals inward and charges 
build up at the center (which is allowed for our compressible 
fluid). 

We note that the fictitious magnetic field bjk/2 
points everywhere in the z direction, its total flux 
through the xy plane being given by 



T = J d(f)dr(rb r< f,) = — J d(j)dr(d^ipd r cos 9) 



= 2tt . 



(83) 

Note that the integrand is just the Jacobian of the map 
from the plane to the sphere defined by the spin-texture 
field: 



(9(v) llp (v)):R 2 



(84) 



This reflects the fact that the fictitious magnetic flux is 
generally a topological invariant, corresponding to the TT2 
homotopy group of the mapping ( 84 1 P3H 



C. Anisotropic resistivity of a 3D spiral 

Consider the texture described by dim — KiZ x m, 
where the spatial rotation stays in the xy plane, but the 
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wave vector n can be in any direction. The spin texture 
forms a transverse helix in the z direction and a planar 
spiral in the x and y directions. Fig. [4] shows such a 
configuration for k pointing along (x + y + z) / v3- The 
fictitious magnetic field b vanishes, but the anisotropic 
resistivity still depends nontrivially on the spin texture: 

7« = [7 + 7?(<9 fc m) 2 ] Sij + • djm 

= (7 + r]K 2 )8ij + ri'KiKj , (85) 

which, according to j = 7~ 1 E, would give a transverse 
current signal for an electric field applied along the Carte- 
sian axes cc, y, or z. 




FIG. 4: (Color online) A set of spin spirals which is topo- 
logically trivial because V# = (and equivalent to the spin 
helix, Fig. [I] up to a global re al-sp ace rotation), hence the 
fictitious magnetic field b, Eq. (76 1, is zero. There is, how- 
ever, an anisotropic texture-dependent resistivity with finite 
off-diagonal components, Eq. (|85l). 



VII. SUMMARY 

We have developed semi-phenomenologically the hy- 
drodynamics of spin and charge currents interacting with 
collective magnetization in metallic ferromagncts, gener- 
alizing the results of Ref. ^2 to three dimensions and 
compressible flows. Our theory reproduces known re- 
sults such as the spin-motive force generated by mag- 
netization dynamics and the dissipative spin torque, al- 
beit from a different viewpoint than previous microscopic 
approaches. Among the several new effects predicted, 
we find both an isotropic and an anisotropic texture- 



classical (topological) Hall effect should be described on 
par with that of the fictitious magnetic field. By calculat- 
ing the dissipation power, we give a lower bound on the 
spin-texture resistivity as required by the second law of 
thermodynamics. We find a more general form, includ- 
ing a term of order (3, of the texture-dependent correction 
to nonlocal Gilbert damping, predicted in Ref. [TT] See 
Eq. ([681). 



Our general theory is contained in the stochastic hy- 
drodynamic equations, Eqs. (58 1, which we treated in 
the highly compressible limit. The most general situ- 
ation is no doubt at least as rich and complicated as 
the classical magnetohydrodynamics. A natural exten- 
sion of this work is the inclusion of heat flows and re- 
lated thermoelectric effects, which we plan to investigate 
in a future work. Although we mainly focused on the 
halfmetallic limit in this paper, our theory is in principle 
a two-component fluid model and allows for the inclu- 
sion of a fully dynamical treatment of spin densities and 
associated flowsP Finally, our hydrodynamic equations 
become amenable to analytic treatments when applied to 
the important problem of spin-current driven dynamics 
of magnetic solitons, topologically stable objects that can 
be described by a small number of collective coordinates, 
which we will also investigate in future work. 
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APPENDIX A: MANY-BODY ACTION 

We can formally start with a many-body action, with 
Stoner instability built in due to short-range repulsion 
between electrons!^ 



S[Vv(r,t),WM)] = j^dtj d 3 r 



(Al) 



dependent resistivity, Eq. ( 53 ) , whose contribution to the 



where time t runs along the Keldysh contour from —00 
to 00 and back. ip a and tp& are mutually independent 
Grassmann variables parametrizing fermionic coherent 
states and ^ + = (^,-0;) and tp = (r/>f, ipi) T ■ The four- 
fermion interaction contribution to the action can be de- 
coupled via Hubbard-Stratonovich transformation, after 
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introducing auxiliary bosonic fields cp and A: 



e iSu/h = exp 



dt / d^rUib^tb^^ 



D[0(r,t),A(r,t)]exp 



eft / d 3 r 



(50 being an arbitrary infinitesimal vector, so that the 
equation of motion is given by m x S m S c s = 0: 

= m x 8 m S eS 



— — tb^tb -\ ■0 + <7"0 



(A2) 



= J^(m x S m a atI ) 



dS 
da„„ 



m x S m F , 



(A5) 



In obtaining this result, we decomposed the interaction 
into charge- and spin-density pieces: 



1 



1 



(A3) 



where Z = J D[ip + ,ip)ei s ^ + ^''f'^ and we have used 
the path-integral representation of the vacuum expecta- 
tion value. a CTAt are the spin-dependent gauge potentials 
Q and F the spin exchange energy, appearing after we 
project spin dynamics on the collective field A. Equa- 



where m is an arbitrary unit vector. It is easy to 
show that (<j>(r,t)) = U(ip + (r,t)ip(r,t)) and (A(r,i)) = 
U(i/j + (r,t)&tjj(r,t)}, when properly averaging over the 
coupled quasiparticle and bosonic fields. 

The next step in developing mean-field theory is to 
treat the Hartree potential 4>(r,t) and Stoner exchange 
A(r, t) = A(r, i)m(r, t) fields in the saddle-point approx- 
imation. Namely, the effective bosonic action 



tion (A5l may be expressed in terms of the hydrody- 



namic variables of the electrons. Defining spin-dependent 
charge and current densities, j£ — {p a ,Sa), by 



Per = 



J (7 = 



dS 
da a 
OS 



= —Re (ip a (-ihV - a a )*Pa) = PaV a , 

(A6) 



5 eff [^(r,t), A(r,t)] = -iHln / D^+^e*^ ^ ;0 ' A) E q- reduces to the Landau-Lifshitz Eq. Min- 



(A4) 

is minimized, SS c s = 0, in order to find the equations 
of motion for the fields <f> and A. In the limit of suf- 
ficiently low electron compressibility and spin suscepti- 
bility, the charge- and spin-density fluctuations are sup- 
pressed, defining mean- field parameters <p and A. Since 
a constant <p only shifts the overall electrochemical po- 
tential, it is physically inconsequential. Our theory is de- 
signed to focus on the remaining soft (Goldstone) modes 
associated with the spin-density director m(r, t), while 
4>(r,t) and A(r, t) are in general allowed to fluctuate 
close to their mean-field values <fi and A, respectively. 
The saddle-point equation of motion for the collective 
spin direction m(r,i) follows from ^miS^rlm] = 0, after 
integrating out electronic degrees of freedom. Because 
of the noncommutative matrix structure of the action 



imizing action (A4| with respect to the <f> and A fields 



gives the anticipated self-consistency relations: 

0(r, t) = U(i> + (r, t)i(r, t)) = U(p + + p_) , 
A(r, t) = U(^+(r, i)a 2 V3(r, *)) = U(p+ -p.). (A7) 

APPENDIX B: THE MONOPOLE GAUGE FIELD 

Let (9, f) be the spherical angles of m, the direction 
of the local spin density, and x<r be the spin up/down 
(er = ±) spinors given by, up to a phase, 



e v sin - 



( A2 ) , it is still a nontrivial problem. The problem sim- 



X-{S, <p) = X+(tt - 9, if + 7r) = 



sin ■ 



(Bl) 



plifies considerably in the limit of large exchange split- 
ting A, where we can project spins on the local magnetic 
direction m. This lays the ground to the formulation dis- 
cussed in Sec. [TTJ where the collective spin-density field 
parametrized by the director m(r, t) interacts with the 
spin-up/down free-electron field. The resulting equations 
of motion constitute the self-consistent dynamic Stoner 
theory of itinerant ferromagnetism. 

In the remainder of this appendix, we explicitly show 
that the semiclassical formalism developed in Sees. [11} 
|IIIB| is equivalent to a proper field-theoretical treatment. 
The equation of motion for the spin texture follows from 
extremizing the effective action with respect to variations 
in m. Because of the constraint on the magnitude of m, 
its variation can be expressed as 5m = 69 x m, with 



The spinors are related to the spin-rotation matrix U(m) 
by Xo =U\<t). The gauge field in m space, which enters 
Eq. (I5JI, is thus given by 

2 \ sin 9 J 

where we used the gradient on a unit sphere: <9 m = 
Ode + y>3 v /sin0. The magnetic field corresponding to 
this vector potential [extended to three dimensions by 
a(m) — > a(#, ip)/m] is given on the unit sphere by 



d m x a" 



r c?e (sin 9a Lf 



ah 



(B3) 
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It follows from Eqs. (|5| and (B2 ) that the spin-dependent 



real-space gauge fields are given by 



-dfj,<p(l — a cos 9) . 



(B4) 



Notice that the a — ± monopole field ( B2 ) , as well as the 



above gauge fields, are singular on the south/north pole 



(corresponding to the Dirac string). This is what allows a 
magnetic field with finite divergence. Any other choice of 
the monopole gauge field (B2 1 would correspond to a dif- 
ferent choice of the spinors ( |B1 1 , translating into a gauge 
transformation of the fields ( |B4| . This is immediately 
seeing by noticing that a™ on (m) — > a™ on (m) + <9 m / CT (m) 



corresponds to a CT „(r, i 



(r,i) + M / ff (m(r,t)). 
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